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^ Abstract 

1—5 

fSJ The impact of right-handed currents in both charged- and neutral-current flavour- violating 

processes is analysed by means of an effective theory approach. More explicitly, we analyse 

,—, the structure of dimension-six operators assuming a left-right symmetric flavour group, com- 

muting with an underlying SU(2)l x SU{2)r x U(1)b-l global symmetry, broken only by two 
ii Yukawa couplings. The model contains a new unitary matrix controlling flavour-mixing in 

the right-handed sector. We determine the structure of this matrix by charged-current data, 
where the tension between inclusive and exclusive determinations of \V u i,\ can be solved. Hav- 
ing determined the size and the flavour structure of right-handed currents, we investigate how 

^-H they would manifest themselves in neutral current processes, including particle-antiparticle 

mixing, Z — > bb, B s ^ — > / u + /i _ , B — > {X s , K, K*}vv : and K — > -nvu decays. The possi- 

:—/ bility to explain a non-standard CP- violating phase in B s mixing in this context, and the 

Q^ comparison with other predictive new-physics frameworks addressing the same problem, is 

also discussed. While a large S^ asymmetry can easily be accommodated, we point out a 

I" — tension in this framework between \V u f,\ and S^k- 

o 

1 Introduction 
> 

^ One of the main properties of the Standard Model (SM) regarding flavour violating processes 



is the left-handed structure of the charged currents that is in accordance with the maximal 
violation of parity observed in low energy processes. Left-handed charged currents encode 
at the level of the Lagrangian the full information about flavour mixing and CP violation 
represented compactly by the CKM matrix. Due to the GIM mechanism this structure 
has automatically profound implications for the pattern of FCNC processes that seems to 
be remarkably in accordance with the present data within theoretical and experimental 
uncertainties, bearing in mind certain anomalies which will be discussed below. 

Yet, the SM is expected to be only the low-energy limit of a more fundamental theory 
in which parity could be a good symmetry implying the existence of right-handed charged 



currents. Prominent examples of such fundamental theories are left-right symmetric models 
on which a rich literature exists. 

Left-right symmetric models were born 35 years ago |Tl-[5j and extensive analyses of many 
observables can be found in the literature (see e.g. pjuland references therein). Renewed 
theoretical interest in models with an underlying SU(2)l x SU(2)r global symmetry has also 
been motivated by Higgsless models 18-11 . However, the recent phenomenological interest 



in making another look at the right-handed currents in general, and not necessarily in the 
context of a given left-right symmetric model, originated in tensions between inclusive and 
exclusive determinations of the elements of the CKM matrix \V u b\ and \V c b\- It could be that 
these tensions are due to the underestimate of theoretical and/or experimental uncertainties. 



Yet, it is a fact, as pointed out and analyzed recently in particular in Ref. [12, 13 , that the 
presence of right-handed currents could either remove or significantly weaken some of these 
tensions, especially in the case of \V u b\- 

Assuming that right-handed currents provide the solution to the problem at hand, there 
is an important question whether the strength of the right-handed currents required for this 
purpose is consistent with other observables and whether it implies new effects somewhere 
else that could be used to test this idea more globally. 

In the present paper we make still another look at the effects of the right-handed cur- 
rents in low energy processes, without specifying the fundamental theory in detail but only 
assuming its global symmetry and the pattern of it breakdown. Specifically we address the 
following questions: 

• What is the allowed structure of the right-handed matrix that governs flavour violating 
processes in the right-handed sector and reduces the tension observed in the inclusive 
and exclusive determinations of |Kj&|? 

• What is the impact of right-handed currents, in combination with SM left-handed 
currents, on particle- antiparticle mixing? 

• Can in this context the known Z — >• bb "anomaly" be solved? 

• Can such effects be seen in rare FCNC decays such as B Sj( i — > /U + /x~ , B — > {X s , K, K*}vi> 
and K — > ttuv! 

As there have been already numerous analyses of right-handed currents in the literature, it 
is mandatory for us to state what is new in our paper: 

• In the spirit of [l4] we use an effective theory approach to describe the effects of right- 
handed currents in flavour violating processes. In fact our work could be considered 
to be a simple generalization of the usual MFV framework to include right-handed 
currents. Indeed the main model-dependent assumption of our analysis is the hypothesis 
that the left-right symmetric flavour group is broken only by two Yukawa couplings. 

• We determine the new unitary matrix V that controls flavour-mixing in the right- 
handed sector by using the data on the tree level charged current transitions u — > d, 
u — > s, b — > u and b — > c and its unitarity. Here the novel feature of our analysis, as 



compared to 12, 15 , is the determination of the full right-handed matrix and not only 



of its b — > c, u elements. 

We point out that the elements of this matrix can be further constrained through the 
FCNC processes, that we analyse in detail. Here our minimalistic assumption about 



the breaking of the flavour symmetry by only two Yukawa couplings plays a key role, 
and distinguishes our work from most of the existing analyses. 

• We point out that the explanation of the different values of \V u b\ from inclusive and 
exclusive semi-leptonic decays with the help of right-handed currents, implying large 
value of \V u b\, strengthens the tension (already existing in the SM) between sin 2/3 and 
SipKs • This tension cannot be solved through the new CP-violating effects in the right- 
handed matrix as the contributions of right-handed currents to B^ — B^ mixing turn 
out to be strongly suppressed due to the ek constraint and the desire to explain the 
enhanced value of the S^ asymmetry in the B s system. 

Our paper is organized as follows. In Section 2 we present the electroweak symmetry 
and the field content of our effective theory. Here we also introduce the right-handed (RH) 
matrix V and discuss some of its properties. In Section 3 we make a closer look at the RH 
charged currents. In Section 4 we use the present knowledge on u — > d, s — > d, b — > c and 
b — >■ u transitions to obtain upper bounds on most of the elements of the matrix V. Section 5 
is devoted to a general discussion of dimension-six operators relevant to neutral currents. 
The impact of these operators is then analysed in particle-antiparticle mixing (Section 6) 
and in processes sensitive to Z- mediated neutral currents (Section 7). A comparison of the 
results obtained in the present framework with those obtained under the MFV hypothesis, 
and in explicit left-right models, is presented in Section 8. We summarize our main results 
and conclude in Section 9. 



2 The model 

2.1 Electroweak symmetry and field content 

The starting point of our analysis is the assumption that the SM is the low-energy limit 
of a more fundamental theory. We don't know the exact structure of this theory, but we 
assume that in the high-energy limit it is left-right symmetric. The difference of left-handed 
(LH) and right-handed sectors observed in the SM is only a low-energy property, due to 
appropriate symmetry-breaking terms. 

In particular, we assume that the theory has a SU(2)l x SU(2)r x U(1)b-l global 
symmetry, explicitly broken only in the Yukawa sector and by the U(l)y gauge coupling. 
Under this symmetry the SM quark fields can be grouped into three sets of LH or RH doublets 
with B — L charge 1/3: 

Q i L = ( u }) , QR=(f) , i = i. ..3. (1) 



\ d L J ' V d R J 

With this assignment the SM hypercharge is given by Y = T^r + (B — L)/2. In order 
to recover the SM electroweak gauge group, we assume that only the SU(2)l and U(l)y 
subgroups of SU(2)l x SU(2)r x U(1)b-l are effectively gauged below the TeV scale. In 
close analogy we can introduce three sets of LH and RH leptons, L l L and L R (including three 
RH neutrinos), with B — L = — 1. 

The breaking of the electroweak symmetry is achieved via the spontaneous breaking 
of SU(2)l x SU(2)r into the vectorial subgroup SU(2)l + r at the electroweak scale. For 



simplicity, we provide an explicit description of this breaking introducing a SM-like Higgs 
field transforming as (2,2) of SU(2)l x SU(2) r , with non-vanishing vacuum expectations 
value (vev): 

H^U L HU R , U l{r) € SU(2) L{R) , (H)=v^ 1 ° 1 ^j. (2) 

Introducing the kinetic term 

^ffi gg s = \Tr[(D,H)^D,H] , (3) 

where D^H = d^H — igW^T a H + ig'HT^B^, we recover the standard tree- level expres- 
sions of W and Z masses for v ~ 246 GeV. However, most of the following discussion on 
flavour-violating effective operators applies as well to models where the SU(2)l x SU(2) r — > 
SU(2)l+r breaking is achieved via a more complicated Higgs sector, or even without a fun- 
damental Higgs field. In the latter case H is replaced by v x U, where U is the unitary field, 
transforming as (2, 2) of SU{2)l x SU(2)r, that encodes the three Goldstone bosons (see 
e.g. Ref. pi). 



2.2 The quark flavour symmetry 

As far as the quark flavour structure is concerned, we assume an SU(3)l x SU(3)r flavour 
symmetry, with the fields transforming as 

Ql(R) -> fL(R)QL(R) , fh(R) € 5 , t/(3) L ( i? ) . (4) 

In order to generate different masses for up- and down-type quarks we introduce two spurion 
fields, P u (d), transforming as UrP u u)Mr under SU(2)r, whose background values break the 
custodial SU(2)l + r symmetry: 

i° Q )< "-(?;)■ (5) 

In addition, we assume that the flavour symmetry is broken by two spurions, both trans- 
forming as (3, 3) under SU(3) L x SU(3) R : 

Y u (d) -> II y u (d) f R ■ (6) 

Finally, we assume an additional U(l) symmetry, under which Y u and Yd have different 
charges, and P u and Pd have the corresponding opposite charges. This symmetry structure 
implies the following invariant quark Yukawa interactions: 

C Y = QlHYuPuQr + Q L HY d PdQ R + h.c. (7) 

This is equivalent to the SM Lagrangian when we take into account the structure of P u and 
Pd- In models with two Higgs doublets the two spurions P u and Pd can be interpreted as the 
remnant of the different vevs of the two Higgs fields. Assuming heavy masses and similar 
vevs for these Higgs fields, this picture cannot be distinguished from the one-Higgs doublet 
case in our effective-theory approach. 



Rotating Ql and Qr in flavour space we can always choose a quark basis where one of 
the two Yukawa couplings is diagonal. We can also rotate the relative phases of the quark 
fields to make this diagonal matrix real. Choosing the basis where Y d is diagonal we can 
write 

i \/2 

Yd | d _ basc = A d , X d = — diag(m d ,m s ,m b ) = diag(y d ,y s ,y b ) , 

*« Id-base = y] ^V, K = diag(m u ,m c ,m t ) =diag(y u ,y c ,y t ) , (8) 

where V and V are two unitary complex 3x3 mixing matrices. V can be identified with the 
CKM matrix, while V is a new unitary mixing matrix that parametrizes the misalignment of 
Y u and Y d in the RH sector. In such a basis, compatible with the SU(3)l x SU(3)r flavour 
symmetry, the mass terms generated by Cy once the Higgs gets a vev are: 

£«— = vu L V^\ u Vu R + vd L \ d d R + h.c. (9) 

In order to diagonalize the mass terms for the up-quarks we need to perform the following 
(flavour-breaking) rotations 

u L ->■ u' L = Vu L , u R ->■ u' R = Vu R , (10) 

with u' L R denoting the mass-eigenstate fields. 

2.3 RH mixing matrix and FCNC spurions 

The new mixing matrix V can be parametrized in terms of 3 real mixing angles and 6 complex 
phases. Adopting the standard CKM phase convention, where the 5 relative phases of the 
quark fields are adjusted to remove 5 complex phases from the CKM matrix, we have no 
more freedom to remove the 6 complex phases from V . In the standard CKM basis V can 
be parametrized as follows 

V = D V %D ] D , (11) 

where Vq is a "CKM-like" mixing matrix, containing only one non-trivial phase and D\j,D 
are diagonal matrices containing the remaining CP-violating phases. For reasons that will 
become clear in the following, it is convenient to attribute the non-trivial phase of V$ to the 
2 — 3 mixing, such that 

C12C13 S12C13 S13 

Vf) =" i -S12C23 - Cl2S23Sl3e~^ C12C23 - S^^S^e' 1 ^ S23Cl 3 e"^ I . (12) 



-C12C23S13 + Sl2S23e 4 ^ -S12C23S13 ~ S23Cl2e^ C23C13 



and 



J 



D v = diag(l,e^,e^3) , D D = diag(e^i, e 1 ^, e**) . (13) 



M) 



Combining the basic spurions Y u and Y d we can build symmetry-breaking terms trans- 
forming as (8, 1) or (1, 8) under the SU(3)l x SU(3)r flavour group. These spurions control 
the strength of FCNCs in the model. Since the only large Yukawa coupling is the top-one, 



only two of such terms are phenomenologically relevant: Y u Yu ~ (8,1) and YuY u ~ (1>8). 
Their explicit expressions in the mass-eigenstate basis of down quarks are 

VirfWi L- b ase = ( V ^l V h « VtV£V 3j , (14) 

(^)^L base = (v^lvh^yt^-^muvohj. (is) 

The YyYu term, which appears in LH mediated FCNCs, has exactly the same structure as 
in the MFV framework [14] . The YdY u term is a new spurion characterizing the strength of 
RH mediated FCNCs in this model. To make contact with the analysis of Ref. |17| , where 
the MFV flavour group (SU(3)q l x SU(3) Ur x SU(3)d R ) with non-minimal breaking terms 
has been considered, our framework corresponds to the introduction of a single non-MFV 
spurion, V, that transforms as (3, 3) under SU(3) UR x SU(3)d R and it is constrained to be a 
unitary matrix. 

3 A first look at the dimension-six operators 
3.1 Preliminaries 

As anticipated, we do not specify the ultraviolet completion of the model. We proceed en- 
coding the effects of the high-energy degrees of freedom by means of the effective Lagrangian 

£ eff = £ S M + ^ J>Of ] , (16) 

where the O- are dimension-six effective operators compatible with the symmetries dis- 
cussed before. Here A is an effective scale, expected to be of 0(1 TeV), and the q are 
dimensionless coefficients. 

In order to build the basis of relevant effective operators, it is first convenient to look at 
the quark bilinear currents invariant under the flavour symmetry defined above. Introducing 
terms with at most two Yukawa spurions (with no more than one Yd), and denoting with T 
a generic Dirac structure, we have 



O(Y ) 
0{Y l ) 
0(Y 2 ) 



QlTQl , QrTQr , (17) 

Q L TY U P U Q R , Q L TY d P d Q R , (18) 

QlTYuY^Ql , QrTY^YuQr . (19) 



Most of these bilinear structures are allowed also in the MFV case. The only two notable 
exceptions are: i) the charged-current component of the RH bilinear QrTQr, and ii) charged- 
and neutral-current components of QrYYuY u Qr. These two will play the central role in our 
paper. We first discuss charged currents, postponing the analysis of neutral currents to 
Section [5j 

3.2 Modification of charged currents 

Our first goal is to understand which operators can probe the rotation in the RH sector, 
namely the flavour-mixing matrix V that appears in the charged-current component of the 



bilinear QrTQr. If we consider operators with only two quark fields, and we ignore RH 
neutrinos (assuming they are heavy), the list of relevant operators is quite small: 

{ Z = Qn-f-nQR L Ll ^L L , 

0f hi = iQr^H^D^HQr , Of = iQirfnQR Tr (h^D^Ht 1 ) . (20) 



Most important, all these operators are equivalent as far as the quark-lepton charged-current 
interactions are concerned. In the case of O r , we generate an effective coupling of the RH 
quark current to the W field after the breaking of the electroweak symmetry: integrating 
out the W leads to a quark-lepton charged-current interaction identical to the one of O r . 
The resulting effective quark-lepton charged-current interaction obtained integrating out 
the W at the tree-level can be written as 

£ eff = ("2M^ + S) ^ dL lhlilVL + F URn " dR lhlilVL + hX ' • (21) 

In the limit cl = cr = we recover the usual SM result. The term proportional to cr is the 
result of the new operators in Eq. (20): cr = —2(cR hl + 2cr H2 — cr (1 ). For completeness, we 



have also included a possible modification of the LH interaction, parametrized by cl- This 



is naturally induced by operators obtained from Eq. ( 20 ) with Qr — > Ql 



In principle, charged-current interactions are potentially sensitive also to operators writ- 



ten in terms of the bilinears in Eqs. (18)-(19). However, as long as we are interested in 



processes where the up- type quarks are of the first two generations, these terms are safely 
negligible, being suppressed by small Yukawa couplings. 



Rotating the up- type fields to the mass-eigenstate basis by means of Eq. (10), and omit- 
ting the prime indices for simplicity, we can finally write 






(1 + c l )VPl + e R VP R d(I L J„VL) + h.c. (22) 



where 



ft-i^a, ^ = 1±^, (23) 

crv 2 v 2 c L v 2 

CR = ~^A2" = A2 (° Rhl Rh2 ~ CRn > ' tL = ~^A^ ' ( ' 

4 Phenomenology of RH charged currents 

In this section we analyse the phenomenology of RH charged currents. In particular, we 
determine the present bounds on the RH mixing matrix V, and we discuss the related impact 



in the determination of the CKM matrix V, using the effective Lagrangian Ctg' in Eq. (22). 
Before starting the phenomenological analysis, we recall that QED and QCD respects 
chiral symmetry. As a result, the two operators in C^q are not mixed by renormalization 
group effects and are multiplicatively renormalized in the same way. This implies that in most 
cases we can incorporate radiative corrections in a straightforward way using SM results. 



4.1 Bounds from u — > d and s — > d transitions 



Within the SM the leading constraints on jV^I are derived from super-allowed (0 + — > + ) 
nuclear beta decays and by the pion decay (ir — > ev) [18]: 



\V ud (0 + - 
\V u d(-ir 



(T 



eu 



icxp 

Ism 

,exp 

Ism 



0.97425(022) 
0.97410(260) 



(25) 
(26) 



By construction, the super-allowed nuclear beta decays are sensitive only to the u — > d 
vector current, while the pion decay is sensitive only the u — > d axial current. As a result, 
the corresponding constraints can be implemented in our effective theory via the conditions 



(l + e L )Vud + e R V ud 
(1 + e L )V ud - e R V ud 



|Kd(0 4 
\V ud (n 



H 



ev 



icxp 

Ism ' 



iexp 

Ism 



(27) 
(28) 



In principle there is also a constraint from the neutron beta decay. However, the situation 
of the neutron life-time measurements is quite confusing at present, and this constraint does 
not add an additional significant new information. 

Since we expect e^/j <C 1, we can expand the above equations to first order in €l,r- 
For simplicity, we also assume €l,r to be real. Solving the above constraints under these 
conditions leads to 



|(1 + ei)K d | = 0.9742 ± 0.0013 , e R Re 



Vud 

V ud 



(0.1 ±1.3) x 10 



-3 



(29) 



The leading constraints on the vector and the axial s — > u currents are derived from 
K — > -kIv and K — >• \xv decays, respectively. Using the SM results obtained in (l9l, 



7T^)| eXP 



SA1 



\V US (K 
\V US (K -)• M^)Ism 

we can impose the following conditions 

1 + e L )V us + e R V u 
[1 + e L )V us - e R V u 

Proceeding as in the u — > d case we find 

|(1 + e L )V us | =0.2248 ±0.0009 , e R Re 



0.2243(12) , 
0.2252(13) , 



\V US {K 



TT^ISM ) 



|Ks(^^^)| oxp 



SM 



Vu 

Vu 



-(2.0 ±3.9) x 10 



-3 



(30) 
(31) 



(32) 
(33) 

(34) 



Since we know that \V, 



u b\ = O(10 3 ), we can neglect |Kj6| 2 m the CKM unitarity relation 
\V ud \ 2 + \V US \ 2 + \V ub \ 2 = 1, and impose that \V ud \ 2 + \V US \ 2 = 1 + O(10~ 4 ). This allows us to 
obtain a determination of e^ at the 10 -3 level starting from the constraints on |(l±£L)V r ud ( s )| 

e L - [(1 + e L ) 2 (\V ud \ 2 + \V US \ 2 + \V ub \ 2 )] V2 - 1 = -(0.2 ± 1.2) x RT 3 . 



inEq. (29) and (34): 



(35) 



Using this result back in Eq. (29) and (34) we finally obtain: 



ud\ 


= 0.9742 ± 0.0013 , 


e R Re(V ud ) = 


(0.1 ±1.3) x 10~ 3 , 


(36) 


us 


= 0.2248 ±0.0009 , 


e R Re(Ks) = 


-(0.5 ±0.9) x 10~ 3 . 


(37) 



Note that if e R = 0(1O -3 ), as expected by naive dimensional analysis for new-physics at the 
TeV scale, the above results do not imply small mixing angles among the first two generations 
in the RH sector. Similar conclusions from the analysis or RH currents in semileptonic K 
decays were reached also in |20|. As we will show in the following, this expectation is 
confirmed also by our analysis of FCNC processes. 



4.2 Bounds from b — > c and b — > u transitions 



Given the smallness of ei derived in Eq. ( 35 ) , we can neglect it in the determination of the 
matrix elements entering b — > c and b — > u transitions, where the best experimental errors 
are at least of 0(1%). 

First we summarize the bounds from b — > c transitions, starting with the consideration of 
the inclusive decay B — > X c lv^. The total rate is easy to handle as the obtained correction, 
being dependent on the fraction of the RH mixing matrix element and the CKM matrix 
element, can be factorized in comparison to the SM tree level decay width. We proceed 
in the following way: the value for V cb , obtained from the comparison of the SM decay 
width and the experimental data, must be equivalent to our effective V cb including the RH 
contributions. We then obtain 



-i/ lincl \ z t/ 

^cfclsM-expj — \ V c 



cb\ 



i + kfll 



Vcb 



'int 



Re e R 



Vcb 

Vcb 



where 121 



l^lsM-ex P = (41-54 ± 0.73) x 10- 
and the strength of the interference rj nt is given by 



Hnt 



m K-) 

" / m b \ 



m b f( 



(38) 



(39) 



(40) 



with f(x) = 1-8x 2 ±8j; 6 



x 



24x 4 log x and h(x) = 1 — 3x 2 + 3x 4 + x 6 + 6(x 2 + x 4 ) log x. 



Numerically r- mt = 0.97 x 10 -3 , thus the impact of the RH current in the inclusive decay 
turns out to be very small. This is consistent with Ref. [15| , where a detailed analysis of the 
inclusive differential distributions has been performed. 

Concerning the exclusive decays B — > D*£ui and B — > D£v£, the consideration of the 
differential decay rate turns out to be more useful as here data has been determined by 
various experiments. Furthermore the consideration of the heavy-quark limit yields an easy 
description, where just one form factor, the Isgur-Wise function, has to be taken into account 
close to the zero-recoil point limit. For completeness, we collect the SM results for the 
differential decay rates close to this kinematical point 



dT(B -> D*tvi) 
dw 



Gp 

48vr 3 



T/ SM I 2 ™ 3 



' cb 



m 



w 



[w 



l) 1/2 P(w)\T(w)f 



(41) 



dT{B 2 m " d = l& 1^1 V* + m D ?mW ~ l? l2 \Q( W )\ 2 , (42) 

where in the -B meson rest frame w = E D (*\/m D (*) and the zero-recoil point limit corresponds 
to w = 1. Here P{w) denotes the phase space factor and J~{w) and Q(w) are the hadronic 
form factors. For w = 1, when the momentum transfer of the leptons is at its maximum, 
-P(l) = 12(ijib — rri£)*) 2 . From a fit of the kinematical distribution around w = 1 the 
experiments determine with high accuracy the products .F(l)|Vd>| and £/(l)|V^,| as well as 
the curvature of the form factors, obtaining 121] , 



^(1)I^*|sm-S, = (35-41 ± 0.52) xl(T 3 , (43) 

<?(l)l^li?_? xp = (42.4 ± 1.6) xlCT 3 . (44) 

In this kinematical limit B — > D*£u£ and B — > Dtvi decays involve only axial and vector 
contributions, respectively. As a result, it is easy to include the RH current contribution. In 
analogy to the inclusive case, hence our conditions read 

l^*li?-S» = |^6-e^c*| , (45) 

\Vcb\sM-e xp = \Vcb + e R V cb \. (46) 

In order to implement these constraints we need to specify the values of the form factors 
at w = 1. Using the lattice determinations Q(\) = 1.074 ± 0.018 ± 0.0016 [22], F(l) = 
0.921 ± 0.013 ± 0.0020 [23], leads to 



\V cb \su-£ P = (39-4 ± 1.7) x 10- 3 , \V cb \^ xp = (38.3 ± 1.2) x 10~ 3 . (47) 



Performing a global fit to V cb and enV c b using the three constraints in Eqs. (39), (45), and 



(46), we then obtain 

lVrf.1 = (40.7 ±0.6) x 10" 3 , e R Re ( -^ ) = (2.5 ± 2.5) x 10" 2 , (48) 

with a modest correlation (p = 0.16). This finally implies 

e R Re(V cb ) = (1.0 ± 1.0) X 10~ 3 . (49) 

In this case the x 2 °f the fit is not good (x 2 /-^dof = 4.3), as also in the SM, because both of 



the exclusive values in Eq. (47) are below the inclusive one. This result cannot be explained 
in terms of RH currents. As pointed out in Ref. [24] , the inconsistency among the different 
determinations of V cb is likely to be due to an overestimate of G(l) on the Lattice. Lowering 
the central value to Q(l) = 0.86, as suggested in [24] , and keeping the same error, leads to 
a much better fit (x 2 /A r dof = 0.9). Since the result for e R Ke(V cb ) obtained in this case is 



perfectly consistent with the one in (49), in the following we will use Eq. (49) as reference 
value. 

We now proceed analysing the constraints from b — > u transitions. As far as the inclusive 
rate is concerned, the structure can be obtained in a straightforward way from the b — > c case 
replacing c — > u. Here the interference term is totally negligible, so we obtain the condition 

(IK^U) 2 = (\V ub \ 2 + \e R \ 2 \V ub \ 2 ) , (50) 



10 



where [21] 

|K^M_exp = (4-11 ± 0.28) x 10- 3 . (51) 

The inclusive determination from B — > nlv, where only the vector current appears, leads to 

|K&lsM-e xp = 1Kb + e R V ub \ = (3.38 ± 0.36) x 1(T 3 , (52) 

where the experimental value is taken from Ref. [21]. Finally, a constraint on the b — > u 
axial current can be obtained from the rare leptonic decay B — > tv. Using the theoretical 
expression 

n(n^^,\ G 2 F m B m 2 T ( m£\ 2 

B(B -)• TU) SM = 1 2" JB\ V ub \ T B , (53) 

8^ V ™> B J 

the experimental result B(B ->• ru) exp = (1.73 ± 0.34) x 10" 4 |25|, and f B = (192.8 ± 
9.9) MeV (26], we get 

\Vub\sM-L P = \Vub~ e R V ub \ = (5.14 ± 0.57) x 1(T 3 . (54) 



As noted first in 12 , here the situation is very favourable for the contribution of RH 
currents, since the axial and vector exclusive determinations are substantially above and 
below the inclusive one (where the interference term is negligible). Performing a global fit 
to V u b and e R V u b using the three constraints we get 

\V ub \ = (4.1 ±0.2) x 10 -3 , e R Re ( -^ ) = -0.19 ±0.07 , (55) 

y V u b J 

with a correlation p = —0.13, namely an evidence of about 2.7cr of a non-vanishing RH 
current contribution. In this case the quality of the fit is excellent (x 2 ~ 0) and substantially 
better than in the absence of RH currents. Most importantly the presence of right-handed 
currents removes the visible discrepancies between the various determinations, as shown in 
Fig. IT] (the SM case corresponds to the top of the vertical axis, where the three determinations 
of \V u b\ are clearly different). 

The above result has been obtained expanding to first order in e R : in this limit we are 
sensitive only to the combination ~R,e(V u b/V u b). Given the non- vanishing result for the RH 
term, we tried also a three-parameter fit, where also Im(V u b/V u b) is free to vary. In this 
case no unambiguous result is found, unless additional conditions are imposed. Imposing the 




condition \V u b\ < \V u b\, the best solution is still the one in Eq. (55), which holds for a large 
interval of Im(V u b/V u b) around zero. Varying the phase of V u b/V u b in a conservative range 
leads to 

|Q?Kfe| = (l-0±0.4)xlO- 3 , for -!I<arg(^| < '- (.",(i) 

4.3 Global fit of the right-handed mixing matrix 

The information we have collected in the previous section can be summarized as follows 

< 1.4 < 1.4 1.0 ±0.4 \ f 3 
\V\~ | - <2.0 I x(-^-) - «•><> 
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Figure 1: Constraints on \V U \\ and e# Re ( y^ ) from d -B — > ir£v (green), B — > X u £u (blue), and 

B — t> tv (orange). The bands denote the ±1<t intervals of the various experimental constraints. 
The ellipse denotes the la region of our best-fit solution. 



where the inequalities correspond to the ±1<t interval and we have assumed small phases 
except for V u b. The entries without figures have very weak direct experimental constraints. 



Altogether the constraints in Eq. (57) seem to be rather weak. However, thanks the unitarity 



condition, they are sufficient to draw a series of interesting conclusions. 

Constraints on e R . The large value of \V u b\ allows us to derive a significant constraint on 
the value of e R from the unitarity of the first row: 



\e R \ = (\e R V ud \ 2 + \e R V us \ 2 + \e R V ub \ 2 Y ' = (1.0 ±0.5) x 10 



(58) 



Given the bound on ei derived in Eq. ( 35 ) , the possibility of ei and e R of the same order 



is perfectly allowed. Note also that the central value in Eq. (58) is in good agreement 



with the naive estimate of models with strong electroweak symmetry breaking, where 
we expect c^jj = 0(1) and A = Anv « 3 TeV. 

We have no information to disentangle the sign of e R and V u b- For simplicity in the 
following we assume e R to be positive. This assumption will not have any consequence 
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for our subsequent analysis since in all observables we always have a similar freedom in 
moving an overall sign from V to the Wilson coefficient of the effective operators. 



Constraints on \Vtb\, \Vt s \, and \Vtd\- Adopting the general parametrization in Eq. (12) 
we find a wide range for the three mixing angles compatible with data. The best-fit 
solution, obtained using only the constraints on e^V^I, e.n|V^ a |, efilKi&l; and e^V^I, 
collected in the previous section, is 






-0.76 


-0.65 


0.88 


-0.31 


0.36 


0.48 


0.57 


-0.67 



F~ 0.88 -0.31 0.36 . (59) 

\ 0.48 0.57 -0.67 / 

However, other solutions with a rather different mixing structure are also compatible 
with data. For instance, the hierarchical scenario with s\s ~ 1 and S12 ~ S23 ~ 0, 
provides also a good fit. What is interesting is that in all cases we can derive non- 
trivial constraints on the elements of the third row, that play a significant role in 
neutral-current observables (to be discussed in the next sections). 
From the unitarity condition on the third column, and the large value of \V u b\, there 
follows a significant constraint on the maximal value of \Vtb\- A large value of \Vtb\ is 
particularly welcome since: i) it minimizes the values of \Vt s \ and \Vtd\, that can induce 
too large contributions to FCNCs; ii) it maximizes the impact of right-handed currents 



in Z — > bb, which could help to improve the agreement with experiments (see Sect. 7.2). 
In the following we thus concentrate on the scenario of maximal \Vtb\- This is achieved 
with the ansatz 

(C12C13 S12C13 S13 

-S12 c 12 I . (60) 

-C12S13 -S12S13 C13 



which corresponds to the C23 — > 1 limit of the general form in Eq. (12). Using the 



ansatz (60), and performing a global fit, we find that at the 90% C.L. 

\Vtb\ < 0.73 , (61) 

where the maximal value is obtained for §13 ~ —0.68 and S12 ~ —0.84. The fact that 
\Vtb\ cannot reach 1 (as in the CKM) implies, in turn, that \Vt s \ and \Vtd\ cannot be 
both vanishing. However, the sensitivity of the fit to S12, which controls their relative 



strength, is very mild. The maximal \Vtb\ in Eq. (61) is obtained for \Vt s \ ~ \Vtd\ ~ 0-5, 
while if we require either \Vts\ or \Vtd\ to be vanishing, then the maximal allowed value 
of \Vtb\ decreases by less than 10%. Only large positive values of §12 are (slightly) 



disfavored by the V us constraint in (37), if 613 is positive 



In summary, we find that the scenario with maximal \Vtb\ is well described by the ansatz 



(|60j) with 

|ci 3 | w -5i3 w 0.7 , e^lx 10" 3 , (62) 

and free S12, provided that sgn(ci3Si2) = — 1. In other words, we find a good description 
with 



m 



±Ci2^ ±S12^ 



C12^ ~S12^ 
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^ ±iJ = -S12 C12 I . (63) 



In the following we assume the matrix Vq as reference structure for the analysis of 
FCNCs. As we will discuss in Section [6j this structure is a key ingredient to generate 
a sizable non-standard contribution to S^. As a result, after we require a large S^ 



(as indicated by recent CDF 27 and DO 28,29] results), most of our conclusions will 



not depend on the choice of this ansatz. It should also be stressed that, contrary to 
the CKM case, having a zero in Vq does not prevent non- vanishing CP-violating effects 



thanks to the extra phases in Eq. (13). 



As far as FCNCs are concerned, the only significant alternative to the ansatz Vq is 
the possibility of a vanishing small \Vtb\, that we can achieve expanding around S13 <C 1 



in (|60j): 

V^ m) = I -S12 ci 2 J . (01) 




As can be seen, this structure is very efficient in escaping all bounds from charged 
currents but for b — > u transitions, where it has a maximal impact. However, it also 
implies naturally small effects in B physics (both in B s and B& meson-anti-meson 
mixing, as well as on rare B decays) and potentially large effects in K physics (which 
are not allowed by data). We thus consider it less interesting with respect to the ansatz 



V (II) in (63). 



4.4 Summary 

We conclude this section with a short summary of the results obtained from the analysis of 
charged currents: 

• RH charged currents can help to reduce the tension between inclusive and exclusive 



determinations in \V u b\, as recently pointed out in 12 . However, contrary to [12] and 
in agreement with [15] , we find that RH charged currents do not have a significant 
impact in the determination of \V c b\. 

The size of the RH charged-current operators necessary to solve the | V u b | problem points 
towards an effective new-physics scale A ~ 3 TeV. 

Thanks to unitarity, the full structure of the RH mixing matrix is quite constrained 
from a global fit of the available constraints. In particular, \Vtb\ is constrained strongly 
by I V u b I and unitarity, which in turn implies non- negligible contributions to FCNCs and 
meson- antimeson mixing through \Vt s \ and \Vtd\- 

The two representative structures of the RH mixing matrix naturally emerge in view 



of the analysis of FCNCs: the structure V in (|63|) is particularly interesting, since it 



could allow large effects in the B s system, without being automatically excluded by the 
tight constraints from the K and Bd systems; the structure Vq in (64) would imply 
sizable effects only in the kaon system. 

5 Dimension-six operators beyond charged currents 

Having constrained the structure of the new mixing matrix from charged current processes, we 
are now ready to analyse the effects generated within our effective theory in neutral-current 
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processes. In particular in the next two sections we analyse the effects in the following set 
of theoretically clean observables: down-type particle-antiparticle mixing (ek and Bd tS -Bd,s 
mixing), rare FCNC decays of B and K mesons with a lepton pair in the final state, and 
Z ->bb. 

The dimension-six effective operators contributing to these processes can be constructed 



by appropriate combinations of the bilinear structures in Eqs. (17)-(19). Rather than pre- 
senting lengthy expressions containing all possible operators, we limit ourself to analyse the 
impact of the most representative ones. In particular, we focus our attention on AF = 2 
operators built in terms of the Q R YuY u r y^Q R bilinear, and AF = 1 operators generating an 
effective right-handed flavour non-universal couplings of the Z boson to down- type quarks. 

5.1 AF = 2 processes 

The complete set of gauge-invariant dimension-six operators contributing to down- type AF = 
2 amplitudes, with the minimum number of Yukawa spurions, is 

0f L = [QUYuY^Qi) 2 , (65) 

0f R = lQ l R (YtYukl»Q J R ] 2 , (66) 

0f R = [QL(Y^)i j ^ I ][Q i R (Y^Y v ) ij ^Q^ , (67) 



The first operator, which is present both in the general MFV framework 14 and in its 
constrained version |30| has been widely analysed in the literature. This operator generates 
short-distance corrections which have exactly the helicity structure and CKM factors of 
the SM short-distance terms. As a result, it cannot modify the SM predictions for the time- 
dependent CP asymmetries in B& and B s decays pO], in particular S^k and S^, respectively. 

A much richer phenomenology is expected from O rr and LR , which are not present in 
the MFV framework and through which the new RH mixing matrix enters the game. For 
this reason, in the following we focus our attention only on these two operators, considering 
the following AF = 2 effective Lagrangian: 

£*™_™0« + «j0». (68) 

Here c RR and clr are flavour-blind dimension-less coefficients, whose size will be discussed 
below. 

5.2 AF = 1 processes 

The list of operators relevant to AF = 1 FCNC processes with a lepton pair in the final state 
can be divided into three categories: 1) operators with two quarks and two Higgs fields; 2) 
operators with two quarks and two lepton fields; 3) dipole-type operators with two quarks 
and one SM gauge field. 

In the first class we have the following two operators 

o%l = iQ R (Y2Y u )^H^D,HQ R , 

0f z2 = iQ i R {YlY u ) in »T i $ R Tr(rfD»HA , (69) 
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and the corresponding LH operators obtained from (69) through Q R — > Ql and Y u •<-)■ Yj 



The latter are of MFV type and have already been analysed in the literature. On the other 



hand, the two operators in (69) give rise to an effective RH coupling of the Z boson of the 
type d R j^ d R Z ^ which is not present in the MFV framework. This coupling is particularly 
interesting since it allows us to establish connections between RH effects in Z — > bb and in 
rare K and B decays. In the following we will analyse these connections by means of the 
effective Lagrangian 

C^ = C -^0f zi + ^Ofl . (70) 

As far as the other two classes of AF = 1 operators are concerned: the operators with 
two quarks and two lepton fields do not lead to effects qualitatively different that those 



obtained after integrating out the Z in Eq. (69). The left-right operators with the photon 



field are tightly constrained by B — > X s j and, similarly to the MFV case, once the B — > Xg'y 
constraint is imposed they do not lead to significant effects in other processes. 



6 Meson anti-meson mixing 



6.1 Preliminaries 

The presence of the right-handed currents in addition to the left-handed ones present in 
the SM can have considerable impact on particle-antiparticle mixing dominantly through 
the generation of the LR operators that renormalize strongly under QCD and in the case 
of K°—K° mixing have chirally enhanced hadronic matrix elements. Below we present the 
general structure of the mixing amplitude M\i from which the observables like ek and the 
CP-asymmetries S^k s an d S^ can be derived. Similarly AM^ and AM S can be calculated. 

6.2 Effective Hamiltonian for AS = 2 transitions 

We will illustrate the full procedure on the example of AS" = 2 transitions. The complete 
set of operators of dimensions six involved in the presence of left-handed and right-handed 
currents consists of the following operators [31] : 



Q\ LL = {SLIM {*L-fd L ) , 

Q\ RR = (SRIM (8jvfd R ) , 

Qi R = (SLl»d L ) (s R rd R ) , 

Qk R = (s R d L ) (s L d R ) , (71) 

where we suppressed colour indices as they are summed up in each factor Pl,r- 

In the SM only Q)[ LL contributes. However, in the presence of right-handed currents 
three additional operators have to be considered. QY renormalizes under QCD, similarly 
to QY LL , fully independently of other operators and in fact the renormalization group running 
for the Wilson coefficient of Q]_ is identical to the one of QY ■ QCD does not care about 
the sign of 75. On the other hand Qf R and Q% R mix under renormalization. 

Let us next assume that at some scale 0(1 TeV), to be denoted by fi R , the new physics 
(NP) is integrated out. In the absence of QCD corrections the effective Hamiltonian for 
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AS = 2 corresponding to fj,R = 0(A) and including only NP contributions reads 

1 
A2 

+Ct R (m,K)Qi R ] . (72) 



[«TMnp = ■h[c 1 ™fa,K)QY LL + cY RR frR,K)QY mt 



The following comments should be made: 

• The coefficients Cf with i = 1,2 and a = VLL,VRR,LR depend generally on the 
system considered (K, B S; d) as is the case also of the SM. 

• At the scale //r, before integrating out the W L fields, the Wilson coefficient C% 
receives only contributions from the NP. Running down to low scales, C\ encodes 
both SM and NP contributions, as discussed below. 

• The coefficient C% (/j,r,K) vanishes at this high scale in the absence of QCD effects 
and it is also vanishing at this scale in the LO renormalization group (RG) analysis. 
However, at NLO it is generally 0(a s (nR)). 



Next usually (72) is evolved by RG down to low energy scales. In this process C^ becomes 
non- vanishing even in the LO approximation. As Q)[ RR and the complex (Qj , 62 ) ^° no ^ 
mix with QY and the RG evolution of the latter operator can be split into SM and NP 
part, the NP contribution to the low energy effective Hamiltonian at a low scale fix can be 
immediately evaluated by means of analytic expressions in [31] . 

While working with Wilson coefficients and operator matrix elements at low energy scales 
is a common procedure, it turns out that for phenomenology it is more useful to work directly 
with Ci(/j,R,K) and with the hadronic matrix elements of the corresponding operators also 



evaluated at this high scale. The latter matrix elements are given by 31 



(K°\Qt\K°) = ^M 2 K F 2 K P?(K), (73) 

where the coefficients P^{K) collect compactly all RG effects from scales below hr as well 
as hadronic matrix elements obtained by lattice methods at low energy scales. 

With this information at hand we can present compactly basic formulae for AS* = 2 and 
AB = 2 observables of interest. 

6.3 Basic formulae for AS = 2 observables 

The off-diagonal element in K° — K° mixing Mf^ can then be decomposed into SM and NP 
parts (RH current contributions in our case) 

MK=(mK) sm +(mE) np , (74) 

where the SM contribution can be found, for instance, in Ref. [31], and the NP part is 
obtained from 

2m K (M^); p = (K°\ [Hf/=ViO] NP |JC°> • (75) 

Using d72b, d73b and d75b we find 



M)np = ^2 m K F% ■ \ (cr LL (^R, k) + cY RR (m, JO) ?? LL (K) 



3A 2 



+CJ: H (m,K)Pt H (K) , (76) 
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parameter 


value 


parameter 


value 


F K 


(155.8 ±1.7) MeV 


26 




AM K 


(5.292 ± 0.009) x lO^ps" 1 (l8) 


Fs d 


(192.8 ±9.9) MeV 


26 




AM d 


(0.507 ± 0.005) ps" 1 36 


Fb s 


(238.8 ± 9.5) MeV 


26 




AM S 


(17.77 ± 0.12) ps- 1 36 


B k 


0.725 ±0.026 : 


26 




\V tb \ 


1±0.06 


37 




B Bd 


1.26 ±0.11 


26 




\Vtd\ 


(8.3 ±0.5) x 10" 3 


37 




Bb s 


1.33 ±0.06 


26 




\Vts\ 


0.040 ± 0.003 


37 




M K 


0.497614 GeV 


18 




sin(2ft) 


0.038 ± 0.003 


37 




M Bd 


5.2795 GeV 


IS 




7 


1.09 ±0.12 


37 




M Bs 


5.3664 GeV 


IS 




exp 
b K 


(2.229 ± o.oiyx 10~ 3 


18 


m t (m t ) 


(163.5 ±1.7) l GeV 


38 




oexp 


0.672 ±0.023 36 





Table 1: Values of the input parameters used in our analysis of AF = 2 processes. 



with 



P{ LL {K) w 0.50, 



p LR 



(K) 



-52 



(77) 



obtained by means of analytic formulae in [31], with the hadronic matrix elements from [32] , 
the updated numerical inputs in Table [TJ and the matching scale [ir = 1.5 TeV. 

The Kl — Ks mass difference and the CP-violating parameter ek are given respectively 



AM K 
£k 



2ReMf 2 + (AM^) LD , 



V2(AM K ) 



(78) 



cxp 



where <p e = (43.51 ± 0.05)° takes into account that tp £ / vr/4 and k £ = 0.94 ± 0.02 (33j[34j 
includes an additional effect from long-distance (LD) contributions (for a recent detailed 
analysis of ek within the SM see [35]). 

Matching the RH effective Lagrangian defined in (68) to the general effective AS = 2 



Hamiltonian in ( 72 ) , the Wilson coefficients at the high scale read 

2 



cY RR (m,K) 



-c RR yfe^-^ 



(VoTtsiVo, 



td 



^rfe*<* 



^W~i 2 ) 2 



-CLRyfe 1 



^-*H*sVu(v y ts (yo) 



td 



cr^(m,K) 



0. 



C -fyte^-^(V t :V td )(c 12 s 12 ), 

(79) 



where the terms on the right-hand side are obtained employing the structure (63) for the RH 
mixing matrix. 

The non-standard contributions to AS = 2 amplitudes are exceedingly large compared 
to the SM term (and compared with data) unless the Wilson coefficients crr and clr or 
one of the two mixing terms c\% or s\% are very small. By construction, crr and clr are 
flavour-blind and therefore the same in the B d and B s system. On the other hand, the c\ 2 
and s\ 2 dependencies in the three systems considered are non-universal, as seen in Table [2j 
with the observables in the if-mixing, B d mixing and .B^-mixing dominated by c\ 2 si 2 , c\ 2 
and S12, respectively. Since both AS = 2 and B d mixing are strongly constrained, and the 



Mixing term 


Matrix 


s ->■ d 


b ->■ rf 


6 — > s 


TO 


CKM 


V t * a V td « -A 5 e"^ 


V^Ftd « A 3 e"^ 


V^VL « -A 2 e"^ 


v t *v tJ 


y(H) 

v m 


|c 12S ~ 12 ^(^-^) 


±|c 12 e*(^-^) 

c 12 g 13 e^3-0?) 





Table 2: Mixing structures relevant to the three down-type AF = 2 and FCNC amplitudes in 
the SM (LH sector) and in the RH sector. In the SM case approximate expressions of the CKM 
factors expanded in powers of A = \V US \ are also shown. In the RH case the two parametrizations 



for the RH matrix are those in Eq. (63) and (64), respectively. 



data from CDF and DO give some hints for sizable NP contributions in the B s mixing, it 
is natural to assume in both scenarios for V that c\i <C 1. In this limit the non- vanishing 



Wilson coefficients relevant for K°-K° mixing and ex at the high scale read with (63) 



cY RR (m,K) 

Ct R (VR,K) 



~yU 2 ^c\ 2 



^^a 



ts v td°l2 



(80) 



where 4> 2 i 



f). Using these expressions in the above formulae for AMk and ek, and 
taking into account the numerical inputs in Table [TJ we obtain: 

(AMx)rh = (AM x ) cxp x [-2.5 x 10 4 x c M 4cos(2^) 

(3 TeV) 2 



far) 



RH 



(AM x ) cxp x ^ 

-1.7 x 10 3 x clrci2 cos(</4 - /3 + p„) 
3.7 xW 6 xc RR c 2 12 sin^; 



\ftA P^ 1 X 



+2.5 x 10 5 x c LR c 12 sin(^ 1 - + /3 S 



A 2 



(3 TeV) 2 



(81) 



(82) 



where, as usual, j3 and /3 S denote the phases of Vm and Vt s in the standard CKM convention: 



V t 



td 



\Vtd\e 



-i/3 



and V ts 



\Vs\e 



-iPs 



(83) 



As can be seen, the RH contribution is potentially very large and, independently of the 
possible value of the CP-violating phase 4> 21 > we § e ^ strong constraints on the combinations 
crrc\ 2 an d clrC\2- In particular, even assuming a vanishing (p 21 , we get 



c RR c\ 2 < 2.0 x 10- 5 , 
clr ci2 < 1.0 x 10~ 6 , 



from (AM K ) R}1 < 0.5{AM K ) cxp , (84) 

from |£k|rh < 0.1|ex| S M • (85) 



While the bound on c RR c^ 2 can only become stronger for non- vanishing values of <j) 2l , the 
constraint on clrC\2 can be relaxed to 3 x 10~ 4 in the fine-tuned scenario where <p 21 cancels 
exactly the CKM phase of V t * s Vtd- As we will show in the next section, these constraints 
imply negligible contribution of RH currents to Bj mixing, which is one of the important 
results of our paper. 
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Before analysing AB = 2 observables, we briefly discuss what happens if do not employ 



the ansatz (63) for the RH mixing matrix and, in particular, if we do not assume £12 <C 1. 



As shown in Table ^1 employing the structure (64) for the RH matrix the mixing structures 
relevant to the kaon system change only by a factor of two. As a result, the corresponding 
bounds on the mixing terms are obtained from (84)-(85) with the replacement cyi —> 2cyiS\2- 



In this case we can escape the kaon bounds and have sizable effects in Bd mixing if £13 is 
not too small. However, the key ingredient for sizable effects in B^ mixing is S12 C 1, a 
configuration that necessarily imply small effects in B s mixing. 

6.4 Basic formulae for AB = 2 observables 

Similarly to the kaon system, for B^ — B^ and B® — B® mixing we can write 



Mh = (M? 2 ) SM + (M? 2 ) Np 



(86) 



where q = d,s. For the NP contribution we find 

1 



« 2 )np 



3A 2 



™ Bq F Bq [ (Cr LL (v R , B) + CY RR (m, B)) Pr L (B) 

+C[ R ( m ,B)P 1 LR (B)V , 



(87) 



with 



P{ LL (B) « 0.70, 



P^ R {B) 



-3.2 



obtained by means of analytic formulae in [31] , with the hadronic matrix elements from [32] , 
the updated numerical inputs in Table [TJ and the matching scale [Ir = 1.5 TeV. 



For the mass differences in the B^ — B® systems we have 



AM q 
Mf 2 



9 I M q I 
z l iw 12l ) 



(Ml 



where 



Mf 2 



12^SM 



SM 



cv 



2i<PB„ 



(q = d,s) 



Mf 2 



SM 



(M 



12JSM 



\{M- 



12/SMl 



„2i£ 



= 2i/3 s 



and Cs q 7^ 1 and ifs q 7^ summarize the NP effects. 



We find then 



and 



Ws 



AM q = (AM q ) SM C Bq 



sin(2f3 + 2ip Bd ) , S^ = sin(2|/3 s | - 2(p Bg 



(89) 

(90) 
(91) 

(92) 
(93) 



with the latter two observables being the coefficients of sin(AM^t) and sin(AM s t) in the time 
dependent asymmetries in B^ — > tpKs and B® — > ijxj), respectively. Thus in the presence of 
non- vanishing ip Bd and ip Bs these two asymmetries do not measure /3 and f3 s but (/3 + y>B d ) 
and (|/3 S | — (fB s ), respectively. 
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The non-vanishing Wilson coefficients at the high scale for the B s d systems are 

2 



C\ (fJR, Bq) 

C-l (fJ'R, B q ) 



-cnnyfe^i-^ 



(YoUVo) 



tq 



-CLRVJe 1 



^^HtVtqiVoUVo)^ ■ 



(94) 



Working in the limit £12 <C 1 (hence S12 ~ 1) we get 

(fjR,Bd)~ ^rVt e 3 c 12 1 L l (fJR,B s ) 

C\ (mh, b s 



qVRRi 
1LR1 



<RR y 4 ( 2i(,<l 



-V2) 



Ct H d>R,B d ) « =F^% 4 e i( ^-^V t * 6 y td c 12 



T 



4 

CL.R 4 



y^e'^s-^v^^ 



V,ty,. 



Ci (vr, B St d) 



0. 



(95) 



where the =F sig n reflects the ± in ( 63 ) . Using these expressions in the formulae for M\ 2 we 
obtain 



(^12)SM+RH = (A^12)SM X 



1 + 



-6.1 x 10 3 x c RR cl 2 e~ 2i ^i+^) ± 4.7 x 10 2 x c LR c l2 e 



i[#§ 1+ 0)\ ( 3 TeV ) S 



A 2 



(M 



12JSM+RH 



(Mf 2 ) SM x 



1 + 



-2.5 x 10 2 x c RR e- 2i ^ + ^ q= 0.9 x 10 2 x cae^^+A) 



(3 TeV) 5 
A^~~ 



(96) 



(97) 



where, similarly to 021; we have defined 0^ = ($5 — <fif). In obtaining the numerical values 
we have evaluated (M^ 2 )sm using the inputs in Table PH 



As anticipated, given the bounds on c RR c\ 2 and cl R c\2 in (84)-(85) following from the 
neutral kaon system, in this framework non-standard contributions to both modulo and phase 
of Bd mixing are safely negligible. On the other hand, sizable contributions to the B s system 
are possible if c RRj lr are in the 10 -3 -10 -2 range and C12 is small enough to satisfy the kaon 
bounds. 



If we do not assume c\i <C 1 and, more generally, go beyond the ansatz (63) for the RH 



mixing matrix, sizable contributions to B& mixing are possible assuming s\i <C 1. However, 
as already stated, this precludes the possibility of large effects in the B s system. 



6.5 Combined fit of Ek and B s mixing 

Here we discuss in more detail the interesting scenario with c rr ^lr m the 10 _3 -10 -2 range 
and C12 <C 1, where we can accommodate a large CP- violating phase in B s mixing, as hinted 



by CDF |27| and DO [28,29 and, at the same time, satisfy the bounds from ek- Values 
of c RRj lr m the 10 _3 -10 -2 range are substantially lower than the 0(1) Wilson coefficients 
determined from charged-current interactions (assuming A = 3 TeV as reference scale). How- 
ever, it is perfectly conceivable that the AF = 2 operators are loop-suppressed with respect 
to the charged-current ones, such that a 10 _3 -10 -2 relative suppression of the corresponding 
Wilson coefficients can naturally be accommodated. 

We stress that a large CP-violating phase in B s mixing is not a clear prediction of the 
model we are considering. As pointed out in the previous section, the only clear prediction 
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is the absence of significant contributions to B^ mixing, implied by the bounds from the 
kaon system, if we require a large CP-violating phase in B s mixing. Still, it is interesting 
to check if this experimental "anomaly" can be solved with reasonable values of the free 
parameters of the effective theory we are considering. To better investigate this point, we 
analyse separately the cases where the leading correction to the SM in B s mixing is induced 
by c RR and c LR , respectively. 

Assuming clr <C crr and neglecting the tiny contribution of (3 S , the constraints of the 
B s system imply 

(AM)sm + rh = _ 2 -**& = (AM.U 

(AM S ) SM (AM S ) SM V ' 

S^ = - ^xl0 2 xc^sin(2^ 2 ) ^ 06±Q3) ^ 

1 - 2.6 x 10 2 x crrc- 21 ^ 

where the latter numerical entry is only a rough indicative value for the CP asymmetry 



favored by the Tevatron experiments (for more details see 27-29 and the model-independent 
analysis in [39] ) . The central values of these equations are full- filled with the following four- 
fold solution: 



CRR S= 


a ±7.3 x 1(T 3 and 


sin(2#j 2 ) p 


a ±0.30 , 


CRR P 


a ±2.3 x 1(T 3 and 


sin(2^ 2 ) p 


a t 0.95 . 



(100) 

As anticipated, these values are in good agreement with the naive expectation of a l/(167r 2 ) ~ 
6 x 10 -3 suppression between this AF = 2 operator (presumably generated at the loop level) 
relative to those contributing to right-handed charged currents (presumably generated at the 
tree level). 

Having fixed crr, we can ask which are the values of the mixing angles necessary to satisfy 
the bounds from the kaon system. To this end, it is first important to note that, due to the 
higher value of \V u b\ determined from charged-current processes with the inclusion of RH 
currents, the prediction of ek in this framework without extra contributions is in excellent 
agreement with data (contrary to what happens in the SM |33|[40|). Indeed the value of 
sin(2/3) determined by tree-level observables only, namely \V u b\ and 7 (following the analysis 
in [37]) is 

sin(2/3)^ e = 0.77 ± 0.05 . (101) 

SM 
tree 



This value is substantially higher than the corresponding result obtained in the SM, sin(2/3) 



0.734 ±0.034 37 , where the inclusive and exclusive determinations of \V u b\ are averaged. As 
a result of this higher value of sin(2/3), the tension between the experimental value of ek and 
its prediction within the SM (see e.g. [35][4l] for a recent analysis) is automatically solved. 

Despite there is no need for non-standard contributions to ek, the theoretical errors on 
this observable allow for extra contributions within ~ ±10% of the SM amplitude. This 
condition is obtained for 

|c 12 ||sin(2^ 1 )| 1/2 < 1.9 x 10~ 3 , for \c RR \ rj 7.3 x 10~ 3 , 

I £12 1 1 sin(2^ 1 )| 1/2 < 3.4 x 10~ 3 , for \c RR \ w 2.3 x 10~ 3 . (102) 

These values are small but not highly fine-tuned: for CP-violating phases of O(0. 1), the 
mixing angle c±2 can reach values of O(10~ 2 ), which are larger than the CKM element \V u b\. 
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Due to the large chiral enhancement of the contribution of Qf R to ek, more fine-tuning 
is required if clr provides the dominant contribution to B s mixing. Indeed repeating the 
above argument for c RR <C clr leads to the following four-fold solution from B s mixing, 

c LR « ±2.0 x 10~ 2 and sin(20f 2 ) « tO.30 , 

(103) 



(104) 



c LR « ±0.6 x 10~ 2 and 


sia(2#j 2 ) « ^0.95 , 


and the following conditions from ek- 




|ci2 tan(<l4i ~ P + &)| < 0-2 x 10~ 4 


for \c LR \ m 2.0 x 10~ 2 


|ci2 sin(^ 1 - p + p„)\ < 0.6 X 10~ 4 


for \c L r\ « 0.6 x 10 -2 



It is clear that in this case the condition on the 1-2 mixing is more stringent than in (102) 



6.6 Summary 

The main results of the AF = 2 analysis can be summarized as follows: 

• A large CP-violating phase in B s mixing, as hinted by the Tevatron experiments, can be 
accommodated for natural values of the free parameters. The two necessary ingredients 
for this mechanism to work are: 1) Wilson coefficients of 0(l/(167r 2 )) for the AF = 2 



operators in (68), assuming as reference scale A = 3 TeV; 2) a RH mixing matrix with 



the following from 



tABs mixing) 




(105) 



where the null entries should not be taken as exact zeros, but rather as very small 
entries. 



According to the RH mixing structure in (105), with improved experimental precision 



it should be possible to resolve the presence of RH currents in s — > u charged-current 
transitions, or a 0(1O~ 3 ) deviation in the determination of \V US \ from K —> iriv and 
K — > lv decays. 

• Thanks to the large value of sin(2/3), following from the inclusion of RH currents in the 
determination on |Vu&|> the prediction of ek in this framework is in excellent agreement 
with data without extra contributions. 

• The combination of a large CP-violating phase in B s mixing and only small NP effects 
allowed by ek implies negligible effects in Bj mixing. This, in turn, implies a tension 
between the measured valued of S^k and the predicted value of sin (2/3) in this frame- 
work. The only possibility to solve this problem is to assume s\2 <^. 1, giving up the 
possibility of NP effects in B s mixing. 

7 Z-mediated FCNCs and Z -> bb 

7.1 Modification of the RH couplings of the Z boson 

As discussed in Section [5j in the AF = 1 sector we focus our attention on the effective 



Lagrangian (70). The effective operators appearing in this Lagrangian are equivalent to 
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3.2, but for the additional insertion of the combination 



On and On analysed in Section 

of Yukawa matrices (Yj_Y u )ij. After the breaking of the electroweak symmetry, they lead to 
the following effective right-handed couplings of the Z boson to down-type quarks: 

9 v 2 {c Rzi +2cr 2 



C 



(Zr) 
eff 



-Vt (VttVtj^R^d^ , 



c w 2A 2 

where cw = cos0w (similarly, in the following we use sw = sinG^y). 

Denoting the effective couplings of the Z to down-type quarks as follows 



C 



eff 



9 

cw 



gi$ I ffJ L + giZ,tf4AZr 



(106) 



(107) 



the SM contribution, evaluated at the one-loop level in the 't Hooft-Feynman gauge in the 
large top-mass limit, is 



(9 j 



zJ SM 



1 1 2 

2 + 3^ 



<% + £^V£V tj C (x t ) , 



nit 



x t 



m 



w 



(9%) 



1 



r)sm 



Sw$ij ■ 



(108) 
(109) 



The loop function Co(xt), that in the large xt limit is gauge independent (g 2 Co(xt) — > 
g 2 xt/8 = y^/4 for m 2 S> m w ), can be found in 42 . Using these notations, the effect of the 



RH operators O r and R can be included as a modification of the RH effective coupling: 



(gt)tot = (g%)sM + (&9%) R H, (Agv 



R )RH 



v2 ( c Rzi + 2c Rz 

2A 2 



-VtVuV tj . 



;no) 



7.2 Z -^ bb 

The experimental determination of the effective couplings of the Z bosons to b quarks result- 
ing from the global fit of electroweak data collected by the LEP and the SLD experiments 

is 



43 



{9L )exp 



-0.4182 ±0.0015 , 
+0.0962 ± 0.0063 . 



;m) 

112) 



(gR Jexp 

While the result for the LH coupling is consistent with the SM prediction, there is a large 
disagreement between data and SM expectation in the RH sector: 

(A<$)ex P = (<$)ex P - (9r)sm = (1.9 ± 0.6) x 10" 2 . (113) 

This deviation could in principle be solved by choosing appropriate couplings for the RH 



operators in (|69|). From the modified RH coupling in (jl 10[) we get 

Jr)rh ? 



(A 5 r 



..eff 
-Zr 



where 



c z R 



(cRzi + 2c Rz 



(3 TeV)^ 



;ii4) 



;ii5) 



and the numerical value has been obtained assuming \Vtb\ 2 

seen, for A = 3 TeV and c Rzi = 0(1), the correction is too small to contribute significantly to 



1/2 (see Sect. 4.3). As can be 



the experimental discrepancy in (113). In principle, the effect could be explained assuming 
A = 1 TeV and cr z1 ,cr Z2 = 0(1). However, as we will show in the following, this possibility 
is ruled out after taking into account the phenomenological bounds from rare B decays. 
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7.3 Rare B and K decays: preliminaries 

General phenomenological analyses about the role of Z-mediated RH currents in B and K 



decays with a lepton pair in the final state can be found in [44 , 45 . While most of the results 
obtained in these two papers can be applied also to the present study, here we go one step 
forward having determined a series of constraints on the flavour structure of the RH mixing 
matrix from other processes. We begin this section by introducing the relevant effective 
Hamiltonians. 

For B S) d — > /i + /x~ channels we generalize the SM effective Hamiltonian to 

4Gf a. 

^eflf = ~ i 2iTg2 V* b V tq X [Y LL (b L ^S L )(fL Llflf , L ) 

+Y L R(b L 'y fJj s L )(fi R >y IM n R ) + Y RL (bRrf8R){fiL r 1 l iV-L) + Y RR (b R ^s R )(p, R 'y^fi R )] , (116) 

where q = d,s. The overall factor allows for an easy comparison with the SM, where 
Yll — Ylr = Yq (xt) and Y RR = Y R l = 0. The contributions to the Y functions with 



RH quark currents obtained by means of the effective Lagrangian ( 106 ) are 



Y RL - Y RR = -T p^ , Y RL + Y RR = -(1- As 2 w ) T Yjjj& , (117) 

where we have defined 



Note that Yll and Ylr are not affected by the effective Lagrangian (106). 

A simple generalization of the SM effective Hamiltonians can be implemented also to 
describe K — > ixvv and B — > {X Sl K,K*jUv decays, with the further simplification that we 
can neglect operators with v R fields, that we assumed to be heavy. In the K — > -kvv case we 
generalize the short-distance effective Hamiltonian to 

n eS = J: 2 V t * s V td x [XLL{K){s L ^d L ) + X RL (K)(s R ^d R )} x {v l1v vl) , (119) 

where the leading SM top-quark contribution yields X R l = and Xll = -^sm = 1-464 ± 
0.041 |46| (for simplicity we omit the sub-leading charm contribution that will be included 
in the phenomenological analysis). With these notations the non-standard contribution ob- 



tained by means of the effective Lagrangian (106) is 



X RL (K) = -T ^ . (120) 

V ts V td 

In the B — > {X s , K, K*}^^ case the general effective Hamiltonian is 

H cS = -^-—^-V? b V ts x [Xll{B s ){Il^sl) + X RL (B s )(b R ^s R )} x (v L >y v v L ) , (121) 
y/2 ^S w 



with 



VIVt, 
X RL (B S ) = -T -*-± . (122) 

V tb V ts 
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7.4 B 



s,d 



fl+fl 



When evaluating the amplitude for B s — > /i + fi by means of (116) the following simplifica- 
tions occur 



(0\h»PR,Ls\B °) = ±-(0\h^ 5 s\B c 



(M^l^ p R,L^\0) = ±-(^1^7^75/210) . (123) 



The resulting branching ratio is then obtained from the known SM expression (see e.g. [47]) 
by making the following replacement 



Y (x t ) ->• Yll + Y RR - Y RL - Y LR = Y t 



tul 



so that 






B{B s ^e + r) = 


--r{B s )^ 

TX 


( a 


U™w 



1117 



F 2 B mjm BsX ll-4^\V t * b V ts \ 2 \Y tot \ 2 



m 



(124) 



(125) 



B, 



The expression for B(B d — > £ + £ ) is obtained by replacing s by d. 

Taking into account that V* b V td » ±ci 2 e^3i/2 and V* b V ts « ±si 2 e i<}!, 32/2 (see Sect . [ZTJ] and 
Sect. 4.3), and using (117), we finally obtain the following expressions for the two branching 



ratios normalized to the SM: 

B(B s ^£ + £-) = B(B s 
B(B d ^£ + £-) = B(B d 



£ + t 



JSM 



JSM 



1 T 7.8 x s 12 e^2 c |ff 
1 ± 37 x ci 2 e 4 ^i c f R 



(126) 



The muon channels are those where the experimental searches are closer to the SM predic- 
tions. The numerical values of the latter, obtained using the relation of B(B q — > jjl + h~~) to 
AMq pointed out in [47], are 



B(B S 



IX' fi 



(3.2 ±0.2) x 10 



-9 



B(B d -> n + pT) = (1.0 ± 0.1) x 10 



-10 



(127) 



These figures should be compared with the 95% C.L. upper limits from CDF [48] and DO [49 
(in parentheses) 



B(B S -»• fi + fj,~) < 3.3 (5.3) x 10" 8 , B(B d -»• fJT(i~) < 1 x 10 



1-8 



Using the results in (|126|) these limits imply 

s 12 cf R <0.54, 



0l20f N 



< 0.30 



(128) 



(129) 



where the bounds have been derived taking into account the interference with the SM (and 
choosing the maximal interference effect). These two limits can be combined to derive the 
following bound 

< 0.62 , (130) 



c z R 



which holds independently of any assumption about the value of ci 2 . Using this bound in 

(Ag^RH < 1 x 10- 3 , (131) 



(114) we get 
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which, by construction, does not rely on any assumption about the value of cyi- It is then 
clear that, despite the presence of a non-standard coupling of the Z boson to RH fermions, 
within our effective theory the constraints from B{B^ S —> £ + £~) prevent a solution to the 
Z — > bb anomaly in (113). 



The bound (130) has been derived from the experimental bounds on both B{B S — > /j + /i~) 
and B{Bd — > /j + /j _ ) in order to show in a simple manner that there is no room for a sizable 
contribution to Z — > bb in our framework, independently of the value of cyi- As far as 
the maximal enhancement of B(B S — > /j + /j~) is concerned, a more stringent bound can be 
derived taking into account the constraints on the effective Z^bn^SR coupling following 
X s £ + t 



from B(B, 



s,d 



44 



45 



In particular, from the bound reported in 45] one finds 



iri x 



v: b vts 



VttVts 



< 1.07 , 



at the 90% C.L. level, or 



^cf R 



< 0.15 . 



(132) 



(133) 



Using this result in (126), the maximal enhancement in B(B S — > /i + /x~) over its SM expec- 
tation does not exceed a factor of 5. This should be contrasted to other NP frameworks, 
in particular to models with non-standard scalar FCNCs, where the present experimental 
upper bound on B(B S — > /_i + /j~) could easily be saturated. 

In case of an 0{\) deviation from the SM in B(B S — > /j + /_i~), a clear prediction of 
our framework, following from the analysis of AF = 2 processes, is the absence of visible 
deviations from the SM in B(Bd — > /i + /_i~). Indeed, as we have seen in Section 6.5, the 
configuration of the RH matrix that could allow to explain a large S^ asymmetry requires 
S12 ~ 1 and C12 < 10~ 2 . When combined with the bound in (133) this condition implies 
negligible non-standard effects in B(Bd — > /i + [i~). 



7.5 B -+ {X s ,K,K*}vi> 

Following the analysis of Ref. [45] , the branching ratios of the B 
the presence of RH currents can be written as follows 



B(B -> Kvv) -- 


= B(B -> Kuu) SM x [1 - 2r/] e 2 , 


B{B ->• K*vv) -- 


= B(B^K*uD) SM x [l + 1.31r/]e 2 


B(B -»■ X s vv) -- 


= B{B -> X s vv)sm x [1 + 0.09?7] e 2 



{X S ,K, K*}vv modes in 



(134) 
(135) 
(136) 



where we have introduced the variables 

2 _ I^llI + |-Xrl| 
e — I Y-SMI2 



-Re(X£ L X RL ) 

|X LL |2 + |X RL |2 ' 



(137) 



in terms of the Wilson coefficient of the effective Hamiltonian ( 12l[ )| 1 | To simplify the nota- 
tions, here and in the following we omit to specify the meson system in the Xlllr functions. 



1 The expressions in Eqs. ( 134 1— ( 136 >, as well as the SM figures in (138), refer only to the short-distance 



contributions to these decays. The latter are obtained from the corresponding total rates subtracting the reducible 



long-distance effects pointed out in 50 
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The updated predictions for the SM branching ratios are (45,50,51 



B(B -> Kvu) m 
B[B -► K*uu) SM 

B(B ->• X s vv) SM 



= (3.64 ± 0.47) x 10 
= (7.2 ±1.1) x 10~ 6 
= (2.7 ±0.2) x 10~ 5 



(138) 



to be compared with the experimental bounds (52T[54 



B{B -► Kvv) 


< 1.4 x 10 


B(B -»■ JfVp) 


< 8.0 x 10 


6(5 -> X s vv) 


< 6.4 x 10 



(139) 



The expressions in Eqs. ( 134 )-( 136 ) are valid for wide class of NP model: all models giving 



rise to the effective Hamiltonian ( 121 ). In our specific framework NP effects are encoded only 



in the RH sector and the Xrl function is given in (122). The variables e and rj then assume 
the following form: 



1 + 



X 2 (x t ) 
T 



>1 



e 2 X (x t ] 



Re 



VttVts 
V t tVts 



1 + 22.1 x | h2cf R | 2 



T 



4.7xai 2 cos(^ 2 )cf a 



,cff 12 



l±22.1x|s 12 4 tt fl 



(140) 
(141) 



Taking into account the bound on si2c|? in Eq. (133), we find that the predictions 
for the exclusive branching ratios can be enhanced by more than a factor of two over the 
corresponding SM estimates. On the contrary, the enhancement in the inclusive mode does 
not exceed 50%. Most important, a clear prediction of RH currents is the anti-correlation of 
the two exclusive modes: if B(B — > Kvv) is enhanced then B(B — > K*vv) is suppressed, and 
viceversajj Not surprisingly, the pattern of these three modes is very similar to what observed 
in the three modes relevant for the determination of \V„h\ in Section 4.2 However, in the 



rare modes the deviations from the SM can in principle be larger than in the charged-current 
decays. 

The correlations among the three B — > {X s , K, K*}vv modes are in principle affected by 
the uncertainty on the CP-violating phase 4>f 2 - However, the same phase enters in S^. If 
we require a large S^ (as hinted by CDF and DO), this uncertainty is strongly reduced, as 
shown in Figure [2} In this plot we show the expectations of the two exclusive branching ratios 



6.5 



for the preferred values of the CP-violating phase (/>g 2 as determined from S^^ in Section 
Two points should be noted: 1) only the modulo of sin(203 2 ) enters in the branching ratios of 
the rare modes (via their cos((/>3 2 ) dependence), as a result, there are only two independent 



choices corresponding to all the solutions in in Eqs. (100) and Eqs. (103); 2) these two 
choices give rise to predictions for the B — > {K,K*}vv branching ratios which are almost 
indistinguishable. As shown in Figure[2], the correlation pattern is very clean and, if observed, 
would provide a clear confirmation of this framework. 



2 In principle, similar correlations could also be established in the B — > K£ + £ and B — > K*£ + £ channels. 
However, in this case the pattern is less clean due to the presence of other effective operators. Moreover, the 



simplifying assumption of considering only the effective Lagrangian ( 70 ) as representative of the dominant RH 
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10 



12 



Br(B^Kvv) x 10~ 



Figure 2: Correlation between B(B — > Kvv) and B(B — > K*vv) in our effective theory. The two 
bands correspond to the two values of | sin(203 2 )| in Eqs. (100), with the uncertainty given by the 
errors on the SM predictions: blue (dark gray) band for |sin(203 2 )| = 0.95, orange (light gray) 
band for | sin(203 2 )| = 0.30. The black point denotes the SM values with the corresponding error 
bars. 



7.6 K — v -kvv 

The SM branching ratios for the two most interesting K 



ttvv modes can be written 



as 



56-59 



B(K A 



TT VV) 



(^) 2 + (^-*-^) : 



where 



-^eff — VtsVtd(XLL + X RL ) 



(142) 
(143) 

(144) 



effects is not necessarily a good approximation for these channels. A detailed analysis of B — > K£ + £ and 
B — > K*£ + £~ in our effective theory goes beyond the purpose of the present paper and we refer to the general 



model-independent analysis in Rcf. 55 . 
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and, as usual, A = \V US \, while k+ = (5.173 ± 0.025) x lQ- n (A/0.225) 8 58 and k l 



(2.29 ± 0.03) x 10 10 (A/0.225) 8 . In the K + case the dimension-six charm quark corrections 



and subleading long-distance effects are characterized by P c = 0.372 ± 0.015 [46,56,59 and 



8Pc,u = 0.04 ± 0.02 57 , respectively. The X e g function can be rewritten as 

X cS = V t * s V td X SM (l + ie ie ) (145) 

where Asm = 1.464 ± 0.041 [46] and we have introduced the two real parameters £ and 6 
that vanish in the SM. The explicit expression for these two parameters in our framework is 

Ze» = ~^ - 5.6 x 10 2 x c l2 ~s 12 e^+e-Mcf R . (146) 

If we ignore the constraints from the AF = 2 processes there is certainly a large room for 



non-standard effects in K — > irvv decays, even taking into account the bound in Eq. (133). 
The situation changes if we implement the constraints on £12 derived from ex, under the 
hypothesis of a large non-standard contribution to S^, analysed in Section 6.5 Here we 
should distinguish two cases: 

• If we can neglect the contribution of Qf to AF = 2 amplitudes, namely if we need 



to implement the bounds in Eqs. (102), then there is room for 0(1) deviations from 
the SM predictions in both rare modes for a sizable range of the phase <ftf 2 . This is 
shown by the red (dark) points in Fig. pi which are obtained imposing this constraint. 
As can be seen, larger deviations from the SM are also possible, but in this case the 
phase (j)f 2 has to be tuned such that sm(2(fif 2 ) ss 0, in order to avoid the ek constraint. 
This fined-tuned configuration gives rise to the two bands of red (dark) points in Fig. [3] 
with large enhancements of only one of the two branching ratios. As noted in [62] this 
structure is characteristic of all NP frameworks where the phase in AS = 2 amplitudes 
is the square of the CP- violating phase in AS" = 1 FCNC amplitudes (this is for instance 
what happens in the Little Higgs model with T parity [63]). 

If the contribution of Qf to AF = 2 amplitudes is dominant, namely if we need to 



implement the stringent bounds in Eqs. (104), then the situation is more constrained 

%2 



Here we can expect a visible deviation from the SM only if the phase <\>\ 2 is tuned such 



that sin(0f 2 — j3 + j3 s ) « 0, where the bounds ( 104 ) become less effective. This give rise 
to the narrow band of green (light) points in Fig. [3] This correlation is very different 
from the one pointed out above, since in this case the leading CP- violating phase in 
AS* = 2 amplitudes is not the square of the CP-violating phase appearing in AS = 1 
FCNC amplitudes. 

If we relax the assumption of sizable NP contributions to S^ the predictions of these two 
modes do not change substantially: the characteristic structures in Fig. [3] are indeed only due 
to the £k constraint. As anticipated, the situation may change only if we could ignore the 
constraints from the AS = 2 processes, assuming the corresponding Wilson coefficients are 
accidentally suppressed. However, we consider this situation highly fine-tuned. Actually it 
should be stressed that the maximal enhancements shown in Fig. [3] also require a considerable 

coefficients. 



amount of fine-tuning, both on the phase 4>f 2 and on the ratio of AS = 2 over AS = 1 Wilson 
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Figure 3: Correlations between B(K + — > tt + vv) and BiK^ — > ifivv) in our effective theory, taking 
into account the ek constraint. The red (dark) and green (light) points are obtained imposing the 
Ek constraint under the assumption of negligible or dominant contribution from Qf , respectively. 

The vertical band correspond to the experimental 



The dashed line is the Grossman-Nir bound 60 
result in [61 



'1 a range) and the black cross to the SM prediction. 



7.7 Summary 



The main results of the AF = 1 analysis can be summarized as follows: 

• The constraints from B s ^ — > /U + /i~ eliminate the possibility of removing the known 
anomaly in the Z —> bb decay with the help of right-handed currents. 

• Contributions from RH currents to B s ^ — > ^ + /i _ , B — > {X s , K, K*}^^, and K — > nvv 
can all be significant, although this is not a general prediction of the model. If the de- 
viations from the SM are sizable, the effects exhibit interesting patterns of correlations. 

• In the B s — > /i + /x~, case the branching ratio can receive an 0(1) enhancement over its 
SM expectation, but it cannot get close to its present experimental bound due to the 
constraint from B — > X s l + l~ . If the RH contribution to S^^ is large, no significant 
enhancement is expected in Bj — > ^ + [i~ . 

• If the RH contribution to S^ is large, the pattern of possible enhancement /suppression 
in B — > {K, K*}uu is unambiguous, as shown in Fig. p] 
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• The pattern of possible enhancement/suppression in the to K — > irvi> modes is largely 
independent of possible RH contributions to Sf</,, but it could help to disentangle the 
case where non-standard AF = 2 amplitudes are dominated by RR or RL operators. 

8 Comparison with MFV and explicit LR models 
8.1 RH currents vs. MFV: general considerations 

The two effective theories are apparently very similar: the low-energy particle content is the 
same and in both cases we have a flavour group broken only by two Yukawas. However, 
the flavour groups are different: SU(3) Ql x SU(3) Ur x SU(S) Dr in the MFV case [II] 
vs. SU(3)l x SU(3)ji in the present framework. The larger flavour symmetry of the MFV 
set-up implies a much more constrained structure. In particular: 

• If the normalization of the two Yukawa couplings is the same as in the SM, as expected 



with a single Higgs doublet, we are in the so-called CMFV regime 30 . In this case 
deviations from the SM are small in all observables. In particular, there is no hope to 
explain a large S^ asymmetry, contrary to what happens with RH currents. 

• The expectation of a vanishingly small S^ in the MFV case remains true also with 
two Higgs doublets and large value of tan /3 if we assume that the Yukawa couplings are 
the only sources of CP violation, as originally assumed in (Ml. In this case the only 
large deviation from the SM is a potentially large B s —> /i + /i~ rate that could easily 
be just below the present experimental bound. As we have seen, in the RH framework 
B s — > /i + /x~ could be enhanced over its SM expectation, but the enhancement cannot 
be as large as in the MFV case at large tan/3 because of the B s — > X s £ + £~ constraint. 

8.2 Comparison with MFV with flavour-blind phases 

The phenomenology of MFV models can be quite different from the case discussed above if 
one relaxes the assumption that the Yukawa couplings are the only sources of CP violation, 
or in the GMFV framework, as denoted in |64| . In particular, it has been recently shown 
that in a two Higgs doublet model (2HDM) with MFV, large tan/3, and flavour-blind CP- 
violating phases, it is possible to generate a large S^ asymmetry and, as a consequence, 



automatically soften the anomalies in Sj,k s and ek in a correlated manner 41 . In this 
set-up the NP contributions responsible for a large S^ are due to the exchange of heavy 
neutral scalars. It is then interesting to compare this solution to the S^ problem with the 
one considered here, where NP contributions originate presumably from new heavy gauge 
bosons. 

Concentrating first on the high energy scales, at which new particles are integrated out 
let us emphasize that the relevant scales of right-handed currents are by roughly an order of 
magnitude larger than the allowed masses of neutral scalars. Also the Lorentz structure of 



the operators generated at the high scale is different in these models. As seen in Eqs. ( 79 ) and 



(95), the presence of RH currents selects from the list in (71 ) non- vanishing initial conditions 
for Qi and Qf , while the leading SU(2)l x U(1)y invariant operator induced by Higgs 
exchanges at the high scale is Q^ R - As a consequence of this different operator structure, 
the NP contributions from RH currents are governed by the RG parameters Pf RR = p^LL 
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and Pi R , while the ones in the 2HDM by P^ R - It is then interesting to observe that model 
independently Pf and P± are roughly of the same magnitude but opposite sign. While 
this sign difference is not relevant in view of unknown signs of the coefficients involved in a 
model with RH currents, it could play some role when a concrete model with RH currents is 
analyzed. 

In the case of the 2HDM with MFV the pattern of NP contributions to K° - K°, B° d - B° d 
and Bg — B® mixings is governed by external quark masses implying naturally largest effects 
in the B® — B®, followed by an order of magnitude 0(md/m s ) smaller effects in B® — B® 
mixing and negligible effects in ek- Still the impact of NP on the B® — B® mixing, uniquely 
following from the requirement of fitting a large S^, softens automatically the ek anomaly 
through the increase of the true value of sin 2/3 resulting from the fit to S^k s ■ It should also 
be noted that the increase of sin 2/3 in the 2HDM automatically favors larger values of \V u b\ 
than those found in B — > irev, but the model does not offer an explanation why the exclusive 
and inclusive determinations give different values of \V U \,\. 

As we have shown in the previous sections, the pattern of deviations from the SM is 
significantly different in the case of RH currents, although some points are common to both 
frameworks. In particular: 

• RH currents provide a natural explanation of the different values of | V u b | following from 
inclusive and exclusive decays, selecting the inclusive determination as the one giving 
the true value of \V u b\. This value is significantly higher than the corresponding value 
obtained from the SM fits and, when combined with the tree level measurement of the 
angle 7, implies sin 2/3 = 0.77 ± 0.05. 

• Similarly to the MFV-2HDM case, the modified determination of sin 2/3 removes the 
ek anomaly even in the absence of direct NP contributions to this observable. 

• In contrast to the MFV-2HDM case, the RH current contributions to B® — B d mixing 
are constrained to be negligible. As a result we should expect S^k s = sin 2/3, implying 
a value of S^k s significantly larger than what determined from experiments. 

• The model has sufficient number of parameters that it can naturally generate a large 
value of S^,^ without any conflict with other data. Contrary to the MFV-2HDM case, 
the large S^ does not imply B s — > £ + £~ close to the present experimental limit. 
Moreover Bj — > £ + £~ receives only small NP contribution so that the MFV relation 
between the the branching ratios for B s — > £ + £~ and B^ — > £ + £~ can be strongly 
violated. 

We conclude then that both models provide interesting solutions to the existing anomalies 
but: 

• the 2HDM with MFV cannot provide the explanation of different values of \V u b\ follow- 
ing from inclusive and exclusive semi-leptonic decays; 

• the model with RH currents provides this explanation but this feature combined with 
tiny contributions to B d — B® mixing implies a ~ 2a discrepancy between the predicted 
value of S^Ks an d its experimental determination. 
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8.3 Comparison with explicit left-right models. 

As anticipated in the introduction, there exist several analyses of flavour observables in 
explicit left-right models (see e.g. J6j[7] and references therein). In all these papers flavour 
observables are a key ingredient to determine the bounds on the masses of the massive 
RH gauge bosons. Most of the existing analyses are focused on the minimal version of the 
model [2}|4], that is characterized by the gauge group SU{2)l x SU(2)r x U(1)b-l, by a 
discrete symmetry connecting the two SU(2) groups, and by the minimal choice for the 
Higgs sector necessary to achieve the two-step breaking SU(2)l x SU(2)r x U(1)b-l — > 
SU(2)l x U(l)y — > C/(l) e .m.- A more general analysis where the latter two hypotheses are 
relaxed can be found in |65|. 

While some features of our effective theory approach can be applied also to these ex- 
plicit models, there are a few important distinctive features of our analysis, related to the 
assumptions behind the Yukawa interactions. 

• In our approach we make no assumptions about the existence of a L-H-R discrete sym- 
metry under which the Yukawa interactions must be invariant. This symmetry, that is 
usually enforced in the minimal models, force the Yukawa matrices to be exactly (or 
approximately) LR symmetric |7|. As a result, V should have a hierarchical structure 
identical (or very similar) to the one of the CKM matrix. While this is a nice feature as 
far as protecting FCNCs, it would prevent us to solve the V u \, problem, which was one 
of the main motivation of our analysis. As far as the structure of V and its impact in 
charged-currents are concerned, our analysis is thus more general than existing analysis 
in explicit LR models. 

• While we do not impose any Lo-R discrete symmetry, we assume a minimal structure for 
the breaking of the flavour group. We assume only two independent Yukawa couplings, 
with an extra protective symmetry forcing them to act only in the up- and down-type 
sector in the dimension- four operators, as shown in Eq. (pi). This extra assumptions have 
been imposed to have an efficient suppression of FCNCs in the left-handed sector, and 
to avoid scalar FCNCs at the tree level. While these assumptions are quite reasonable 
and can easily be implemented in explicit LR models, they do not represent the most 
general possibility. It is worth to stress that some of our phenomenological conclusions, 
such as the absence of NP effects in B^ mixing, after we require large NP effects in B s 
mixing, do depend on this assumption. 

9 Conclusions 

The possibility that at very short distance scales the nature is left-right symmetric appears 
to be intriguing. As at low energy scales the parity is maximally broken and charged weak 
interactions exhibit left-handed structure, the right-handed weak currents, if present in na- 
ture, must be coupled to new heavy gauge bosons that are at least by one order of magnitude 
heavier than the W^ and Z bosons of the SM. 

While such heavy gauge bosons could be discovered at the LHC in the coming years, they 
can also manifest themselves in low energy processes. In our model they are represented 
by new effective operators containing right-handed currents with their Wilson coefficients 
encoding the information about the fundamental theory, in particular the relevant couplings. 
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In the present paper we have analyzed the impact of right-handed currents in both 
charged- and neutral-current flavour-violating processes by means of an effective theory ap- 
proach. To this end we have assumed a left-right symmetric flavour group, commuting with 
an underlying SU{2)i x SU(2)r x U{1)b-l global symmetry, broken only by two Yukawa 
couplings. Having identified the leading six dimension operators in this model, we performed 
a rather detailed analysis of those low energy observables, which could help to support or 
falsify the presence of this NP at scales probed soon directly by the LHC. 

The central role in our model is played by a new unitary matrix V that controls flavour- 
mixing in the right-handed sector. Using the data on the tree level charged current transitions 
u ^ d, n— ^s,6— ^n and b — > c and the unitarity of V we could determine the structure of 
this matrix and demonstrate, following other authors, that the tension between inclusive and 
exclusive determinations of \V u b\ can be solved with the help of right-handed currents. The 
resulting true value of \V u b\ turns out to be (4.1 ± 0.2) x 1CP 3 , in the ball park of inclusive 
determinations. The novel feature of our analysis as compared with [12] is the determination 
of the full right-handed matrix, and not only its selected elements, making use of unitarity. 
We also find that while RH currents are very welcome to solve the u \V u b\ problem" they do 
not have a significant impact on the determination of \V c b\ (as also pointed out in |15|). 

Having determined the size and the flavour structure of right-handed currents that is con- 
sistent with the present data on tree level processes and which removes the "|V^&| -problem", 
we have investigated how this NP would manifest itself in neutral current processes, includ- 
ing particle-antiparticle mixing, Z — > bb, B S( i — > |U + //~, B — > {X s ,K,K*}vv and K — > ttuu 
decays. We have also addressed the possibility to explain a non-standard CP-violating phase 
in B s mixing in this context, and made the comparison with other predictive new-physics 
frameworks addressing the same problem. 

The main messages from this analysis are as follows: 

• The presence of RH currents in the model, in conjunction with the already present SM 
left-handed currents generates in addition to the operators with (V + A) x (U + A) Dirac 
structure, also left-right operators (V — A)x(V + A). The contributions of the latter are 
known to be strongly enhanced at low energies through renormalization group effects 
and in the case of Ek and AMk through chirally enhanced hadronic matrix elements of 
(V — A) x (V + A) operators. Consequently these observables put severe constraints on 
the model parameters ( as also known from various studies in explicit LR models |6|). 

• The desire to generate large CP-violating effects in _B s -mixing, hinted for by the en- 
hanced value of S^ observed by the CDF and DO collaborations, in conjunction with 
the e^-constraint, implies additional constraints on the shape of V. In particular 
C12 *C 1 and consequently §12 ~ 1. The pattern of deviations from the SM in this 
model is then as follows. 

• The S^ and Ek anomalies can be understood. 

• As a consequence of the large value of su, it should be possible to resolve the presence 
of RH currents also in s — > u charged-current transitions. Here RH currents imply a 
0(1O -3 ) deviation in the determination of \V US \ from K — > tt£v and K — > £u decays. 

• The "true value" of sin 2/3 determined in our framework, namely the determination of 
the CKM phase ft on the basis of the tree-level processes only, and in particular of 
\V u b\, is sin 2/3 = 0.77 ± 0.05. This result is roughly 2a larger than the measured value 
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S-tpKs = 0.672 ± 0.023. This is a property of any explanation of the u \V u b\ -problem" 
by means of RH currents, unless the value of \V u b\ from inclusive decays will turn out 
to be much lower than determined presently. In general, such discrepancy could be 
solved by a negative new CP-violating phase in B^ — B^ mixing. However, we have 
demonstrated that this is not possible in the present framework once the ek constraint 
is imposed and large S^ is required. Thus we point out that simultaneous explanation 
of the "|T4b| -problem" and of S^k s = 0.672 ±0.023 is problematic through RH currents 
alone. 

• The present constraints from B s ^ — > n + n~ eliminate the possibility of removing the 
known anomaly in the Z — > bb decay with the help of right-handed currents. On top 
of it, the constraint from B — > X s l + l~ precludes B s — > [J> + fJ>~ to be close to its present 
experimental bound. Moreover NP effects in B^ — > £ + £~ are found generally smaller 
than in B s -> l + T . 

• Contributions from RH currents to B — > \X Sl K,K*}vv and K — > isvv decays can 
still be significant. Most important, the deviations from the SM in these decays would 
exhibit a well-defined pattern of correlations. 

We have compared this NP scenario with the general MFV framework and with more 
explicit NP models. Particularly interesting is the comparison with the 2HDM with MFV, 
large tan/3, and flavour-blind CP-violating phases, where the S^s and ek anomalies can also 



be accommodated 41 . What clearly distinguishes these two models at low-energies is how 
they face the "\V u b\ -problem" (which can be solved only in the RH case) and the "sin 2/3-S^ k 
tension" (which can be softened only in the 2HDM case) . But also the future results on rare 
B and K decays listed above could in principle help to distinguish these two general NP 
frameworks. 

Restricting the discussion to these two NP frameworks, it appears that a model with an 
extended scalar sector and right-handed currents could provide solutions to all the existing 
tensions in flavour physics simultaneously. This possibility can certainly be realized in explicit 
left-right symmetric models, where an extended Higgs sector is also required to break the 
extended gauge symmetry. However, these extensions contain many free parameters and 
clear cut conclusions on the pattern of flavour violation cannot be as easily reached as it was 
possible in the simple framework considered here and in |41| . 

We are looking forward to the upcoming experiments at the LHC, future B factories, and 
rare K decay experiments, that should be able to shed more light on the role of right-handed 
currents in flavour physics. 

Acknowledgments 

We thank Wolfgang Altmannshofer and Tillmann Heidsieck for useful discussions. G.I. would 
like to thank Jernej Kamenik and Martin Gonzalez- Alonso for useful discussions in the early 
stage of this work. We all thank the Galileo Galilei Institute for Theoretical Physics for 
the hospitality and partial support during the completion of this work. This research was 
partially supported by the Cluster of Excellence 'Origin and Structure of the Universe', by 
the Graduiertenkolleg GRK 1054 of DFG, by the German 'Bundesministerium fiir Bildung 
und Forschung' under contract 05H09WOE, and by the EU Marie Curie Research Training 
Network contract MTRN-CT-2006-035482 (Flavianet). 



36 



References 

[1] J. C. Pati and A. Salam Phys. Rev. Lett. 32 (1974) 1083. 

[2] R. N. Mohapatra and J. C. Pati Phys. Rev. Dll (1975) 2558. 

[3] R. N. Mohapatra and J. C. Pati Phys. Rev. Dll (1975) 566-571. 

[4] G. Senjanovic and R. N. Mohapatra Phys. Rev. D12 (1975) 1502. 

[5] G. Senjanovic Nucl. Phys. B153 (1979) 334. 

[6] Y. Zhang, H. An , X. Ji, and R. N. Mohapatra Nucl Phys. B802 (2008) 247-279, 
jarXiv: 0712 .4218] . 



[7] A. Maiezza, M. Nemevsek, F. Nesti, and G. Senjanovic |arXiv : 1005 . 5160 

C. Csaki, C. Grojean, L. Pilo, and J. Terning Phys. Rev. Lett. 92 (2004) 101802, 
|hep-ph/0308038|. 

[9] Y. Nomura JEEP 11 (2003) 050, [ hep-ph/0309189| . 

[10] R. Barbieri, A. Pomarol, and R. Rattazzi Phys. Lett. B591 (2004) 141-149, 



|hep-ph/0310285 



[11] H. Georgi Phys. Rev. D71 (2005) 015016, |hep-ph/0 408067|. 



[12] A. Crivellin Phys. Rev. D81 (2010) 031301, arXiv: 0907. 2461 



[13] C.-H. Chen and S.-h. Nam Phys. Lett. B666 (2008) 462-466, [arXiv:0807.0896| . 

[14] G. D'Am brosio, G. F. Giudice, G. Isidori, and A. Strumia Nucl. Phys. B645 (2002) 
155-187, [hep-ph/0207036]. 



[15] R. Feger, V. Klose, H. Lacker, T. Lueck, and T. Mannel arXiv: 1003.4022 



[16] G. Isidori |arXiv : 0911 ■ 3 219 . 



[17] T. Feldmann and T. Mannel JEEP 02 (2007) 067, |hep-ph/06fl 095]. 

[18] Particle Data Group Collaboration, C. Amsler et. al. Phys. Lett. B667 (2008) 1. 



[19] M. Antonelli et. al. |arXiv: 1005.2323[ 



[20] V. Bernard, M. Oertel, E. Passemar, and J. Stern JEEP 01 (2008) 015, 
arXiv: 0707. 4194]. 



[21] M. Antonelli et. al. arXiv: 0907. 5386 



[22] M. Okamoto PoS LAT2005 (2006) 013, [hep-lat/0510113 



[23] C. Bernard et. al. Phys. Rev. D79 (2009) 014506, 



arXiv: 0808. 2519 



[24] P. Gambino, T. Mannel, and N. Uraltsev arXiv: 1004.2859. 



[25] UTfit Collaboration, M. Bona et. al. Phys. Lett. B687 (2010) 61-69, 
|arXiv: 908. 3470], 

[26] J. Laiho, E. Lungh i, and R. S. Van de Water Phys. Rev. D81 (2010) 034503, 
arX i~v:0910.2928"] . 

[27] CDF Collaboration, T. Aaltonen et. al. Phys. Rev. Lett. 100 (2008) 161802, 
|arXiv:0712.2397] . 



37 



[28] DO Collaboration, V. M. Abazov et. al. Phys. Rev. Lett. 101 (2008) 241801, 
I arXiv: 0802. 22551. 



[29] DO Collaboration, V. M. Abazov et. al. arXiv: 1005.2757. 



[30] A. J. Buras, P. Gambino, M. Gorbahn, S. Jager, and L. Silvestrini Phys. Lett. B500 
(2001) 161-167, |he p-ph/00070 85|. 

[31] A. J. Bu ras, S. Ja ger, and J. Urban Nucl. Phys. B605 (2001) 600-624, 
|he p-ph/010 2316|. 

[32] R. Babich et. al. Phys. Rev. D74 (2006) 073009, [hep-lat/0605016]. 

[33] A. J. Buras and D. Guadagnoli Phys. Rev. D78 (2008) 033005, [arXiv: 0805 .3887] . 

[34] A. J. Buras, D. Guadagnoli, and G. Isidori 



arXiv: 1002.3612 



[35] J. Brod and M. Gorbahn arXiv: 1007.0684 



[36] Heavy Flavor Averaging Group Collaboration, E. Barberio et. al. arXiv: 0808. 1297 



[37] UTfit Collaboration, M. Bona et. al. JEEP 03 (2008) 049, [arXiv: 0707. 0636 



[38] Tevatron Electroweak Working Group Collaboration arXiv : 0903 . 2503 



[39] Z. Ligeti, M. Papucci, G. Perez, and J. Zupan arXiv: 1006.0432 



[40] E. Lunghi and A. Soni Phys. Lett. B666 (2008) 162-165, |arXiv : 0803^ 4340 



[41] A. J. Buras, M. V. Carlucci, S. Gori, and G. Isidori |arXiv : 1005 . 5310 

[42] G. Buchalla, A. J. Buras, and M. E. Lautenbacher Rev. Mod. Phys. 68 (1996) 
1125-1144, |nep-ph/9512380|. 

[43] LEP and SLD Electroweak Working Groups Collaboration Phys. Rept. 427 (2006) 



257, [hep-ex/0509008 



[44] G. Buchalla, G. Hiller, and G. Isidori Phys. Rev. D63 (2000) 014015, 
|he p~ph/00061 36|. 

[45] W. Altmannshofer, A. J. Buras, D. M. Straub, and M. Wick JEEP 04 (2009) 022, 
|arXiv: 902 .0160]. 

[46] A. J. Buras, M. Gorbahn, U. Haisch, and U. Nierste JEEP 11 (2006) 002, 
|he p-ph/060 3079|. 

[47] A. J. Buras Phys. Lett. B566 (2003) 115-119, |h ep-ph/0303060| . 

[48] CDF Collaboration, T. Aaltonen et. al. Phys. Rev. Lett. 100 (2008) 101802, 
ar XTv: 0712. 1708] . 

[49] D O Collabora tion, V. M. Abazov et. al. Phys. Rev. D76 (2007) 092001, 
|arXiv : 0707. 3997 . 



[50] J. F. Kamenik and C. Smith Phys. Lett. B680 (2009) 471-475, [arXiv: 0908. 1174 



[51] M. Bartsch, M. Beylich, G. Buchalla, and D. N. Gao JEEP 11 (2009) 011, 



arXiv: 0909. 1512 . 



[52] ALEPH Collaboration, R. Barate et. al. Eur. Phys. J. C19 (2001) 213-227, 
|hep-ex/0010022|. 



38 



[53' 

[54 

[55 
[56; 

[sr; 

[58; 
[59; 

[60 

[61 

[62 
[63 
[64 



BELLE Collaboration, K. F. Chen et. al. Phys. Rev. Lett. 99 (2007) 221802, 
[arXiv:070T .0138|. 

BABAR Collaboration, B. Aubert et. al. Phys. Rev. D78 (2008) 072007, 
arXiv: 0808. 1338]. 



W. Altmannshofer et. al. J 'HEP 01 (2009) 019, [a rXiv: 0811. 1214]. 

A. J. Buras, M. Gorbahn, U. Haisch, and U. Nierste Phys. Rev. Lett. 95 (2005) 
261805, |hep-ph/0508165|. 

G. Isidori, F. Mescia, and C. Smith Nucl. Phys. B718 (2005) 319-338, 



|hep-ph/0503107 



F. Mescia and C. Smith Phys. Rev. D76 (2007) 034017, arXiv : 0705 . 2025| 



J. Brod and M. Gorbahn Phys. Rev. D78 (2008) 034006, arXiv:0805.4119 



Y. Grossman and Y. Nir Phys. Lett. B398 (1997) 163-168, [hep-ph/9701313 



E949 Collaboration, A. V. Artamonov et. al. Phys. Rev. Lett. 101 (2008) 191802, 
arXiv: 0808. 2459]. 



M. Blanke Acta Phys. Polon. B41 (2010) 127, arXiv: 0904. 2528 



M. Blanke et. al. JEEP 01 (2007) 066, [hep-ph/0610298 



A. L. Kagan, G. Perez, T. Volansky, and J. Zupan Phys. Rev. D80 (2009) 076002, 



arXiv: 0903. 1794]. 
[65] P. Langacker and S. Uma Sankar Phys. Rev. D40 (1989) 1569-1585. 



39 



